GSI-2004-13, IPNO-DR-04-03 



High-resolution velocity measurements on fully identified light nuclides 
produced in ^^Fe + hydrogen and ^^Fe + titanium systems 
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New experimental results on the kinematics and the residue production are obtained for the 
interactions of ^ Fe projectiles with protons and "''^Ti target nuclei, respectively, at the incident 
energy of 1 yl GeV. The titanium-induced reaction serves as a reference case for multifragmentation. 
Already in the proton-induced reaction, the characteristics of the isotopic cross sections and the 
shapes of the velocity spectra of light residues indicate that high thermal energy is deposited in the 
system during the collision. In the ^^Fe+p system the high excitation seems to favour the onset of 
fast break-up decays dominated by very asymmetric partitions of the disassembling system. This 
configuration leads to the simultaneous formation of one or more light fragments together with one 
heavy residue. 

PACS numbers: 25.40. Sc, 25.70.Pq, 24.10.-i, 21.10.Gv 



INTRODUCTION 



For the last decades, the investigation of the maxi- 
mum excitation energy that a nuclear system can hold 
has remained as much a challenge as the description of 
the decay of a hot collision remnant, excited beyond the 
limits of nuclear binding. It is commonly assumed that 
other decay modes than fission and evaporation prevail 
at high excitation energy. These modes are often de- 
scribed as a simultaneous break-up of the hot system 
in many parts, named "multifragmentation". The exci- 
tation energy above which multifragmentation appears 
is still a source of intense theoretical and experimental 
research. A point of particular interest is to recognise 
the distinguishing traits denoting this decay mode when 
the excitation is just sufficient for its onset. In line with 
this investigation, one foremost aspect of intense discus- 
sion is the connection of the kinematics of the residues 
to the kind of equilibration process involved in the ear- 
liest stages of the decay. This question is related to the 
complementary effort in constructing physical models to 
deduce the formation cross sections of the residues when 
the excitation energy of the system is taken as initial 
condition. Especially light residues are suited for this 
purpose. Several details of the deexcitation mechanism 
could emerge from the kinematics of fight fragments, due 
to the high sensitivity in probing the Coulomb field of the 
decaying system. Moreover, the distribution of their iso- 
topic cross sections carry additional signatures connected 
to different decay modes. 
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The experimental data exploited in this work are part 
of a vast experimental campaign devoted to the collec- 
tion of nuclear data for the design of accelerator-driven 
subcritical reactors |l| and to the investigation of spalla- 
tion reactions in the cosmos J2| • The full set of produc- 
tion cross sections and fragment mean recofi velocities 
measured in this framework for the reaction ^^Fe-|-p and 
56pg_^natrp- ^^^ energies ranging from 300 to 1500 MeV is 
presented elsewhere J3|. 



A. The formation of light residues 

Light residues can be generated in several kinds of pro- 
cesses. One of these, the binary decay of an excited 
greatly thermalised complex, named compound nucleus, 
was widely studied |j|. We might also recall that evap- 
oration of nucleons and light nuclei and symmetric fis- 
sion are just the opposite extremes of the manifestation 
of this process: there is a gradual transition from very 
asymmetric to symmetric configurations in the division 
of decaying compound nuclei, and thus all binary decays 
of a greatly thermalised system can be named fission in 
a generalized sense. This generalization was introduced 
by Moretto @, la|- A compound system far below the 
Businaro-Gallone point p, |^ (like iron-like nuclei) un- 
dergoes very asymmetric fission, resulting in a character- 
istic U-shape in the mass distribution of the yields. A 
minimum located at symmetry in the yield mass spec- 
trum corresponds to a maximum placed at symmetry in 
the ridge lines of the potential. In configurations where a 
heavy partner is present, the whole decay process is domi- 
nated by the binary decay, and an additional evaporation 
of single nucleons would not disturb the kinematics re- 
markably. Such a process exhibits the typical feature of 
the population of the shell of a sphere in velocity space, 
in the reference frame of the mother nucleus. 



At high excitation, multifragmentation becomes the 
competing process to compound-nucleus reactions. 
There is a fundamental difference between the binary de- 
cay of a compound nucleus and the simultaneous disinte- 
gration of a hot collision remnant in several constituents. 
The difference is in the kind of instabilities which are the 
reason for the decay, and is reflected in the kind and in 
the time evolution of the consequent equilibration process 
followed by the system. 

A hot nucleus with an excitation energy above the 
threshold for emission of particles or clusters (includ- 
ing fission) has the possibility to decay by any of the 
open channels. If the excited system is not too hot, the 
favoured process is a reordering of its configurations: a 
great number of arrangements are available where all nu- 
cleons remain in states below the continuum, occupying 
excited single-particle levels around the Fermi surface. 
Oscillations in fission direction are included in this pic- 
ture as well, but too rarely the fission barrier is reached. 
Rather seldom, compared with this thermal chaotic mo- 
tion of the system, one nucleon acquires enough energy 
to pass above the continuum and may eventually leave 
the nucleus. This picture might be extended to cluster 
decay and to fission. Since this decay is a rare process, 
one evaporation event, or fission event, proceeds after the 
other, sequentially. In this process, the compound sys- 
tem follows a dynamic trajectory in deformation space, 
which is governed by the potential-energy surface and the 
dynamic properties of the compound system, related, for 
instance, to the inertia tensor and dissipation tensor. All 
decays are binary. 
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FIG. 1: Idealistic plot of the phase diagram of nuclear mat- 
ter, deduced from a Skyrme force f |27l | parameterised accord- 
ing to 123). Pressure is shown as a function of the average 
relative nucleon distance r normalized to the distance ro at 
ground state. System configurations are drawn as possible 
final results of the expansion phase. When the thermaliza- 
tion path leads to the coexistence region, out of the spinodal 
region, damped density fiuctuations occur. In the spinodal 
region density fiuctuations are unstable and lead to cracking. 
At low density freeze-out is attained with different possible 
partition configurations: fragments are free to leave the sys- 
tem. 



If the system becomes drastically more unstable, this 
picture is not vaHd anymore. The exploration of possible 
states of the excited system includes numerous unsta- 
ble configurations. Thus, the disintegration can not be 
understood as a sequence of binary decays, but rather 
portrayed as a simultaneous break-up in several con- 
stituents ja, llO, ini U2> UM- The disintegration is si- 
multaneous in the sense that it evolves in so short a time 
interval (lO^'^^-lO^^^s) that the ejected fragments can 
still exchange mutual interactions during their accelera- 
tion in the Coulomb field of the system. In heavy-ion 
collisions, part of the excitation could be introduced in 
the system in the form of compressional energy. Accord- 
ing to the impact parameter and the incident energy, the 
interaction might result in a very complex interplay be- 
tween dynamic effects (beside compression, also defor- 
mation and rotation degrees of freedom) and thermal 
excitation. This is the case of central collisions in the 
Fermi-energy range. On the contrary, peripheral heavy- 
ion collisions at relativistic energies may be rather pic- 
tured according to an "abrasion" process [ij, [T^ , where 
the remnant is formed by the spectator nucleons, heated 
by mainly thermal energy. In this case, the role of com- 
pressional energy has minor incidence. Even with proton 
projectiles, the multifragmentation regime might be ac- 
cessible when very high excitation is introduced in the 
nucleus. In reactions induced by relativistic protons (but 
also by very light nuclei), the dynamic effects of the 
collision have even smaller importance. The excitation 
energy is almost purely thermal. Some authors even 
attributed the specific name of "thermal multifragmen- 
tation" to this particular process (see the review arti- 
cles [lall3)- It might be suggested that proton-induced 
relativistic coUisions are better suited than ion-ion colli- 
sions for investigating thermal properties of nuclear mat- 
ter (e.g. [H IH lia El Ea)- In finite nuclei, the 
transition from the fission-evaporation mode to multi- 
fragmentation manifests rather smoothly. This open- 
ing of break-up channels even inspired interpretations 
in line with the liquid-gas phase transition of nuclear 
matter IHIl^, El Hi ■ The similarity of the nucleon- 
nucleon interaction with the Lennard-Jones molecular 
potential suggests that infinite neutral nuclear matter re- 
sembles a Van-der-Waals fiuid J23|. As shown in fig. Q] 
also in the phase diagram of nuclear matter an area of 
liquid-gas coexistence can be defined. In this region, the 
"dense" phase of nuclear droplets is in equilibrium with 
the "gaseous" phase of free nucleons and Hght complex 
particles. Within the Hartree-Fock approximation, ac- 
cording to the type of Skyrme force chosen for obtaining 
the nuclear equation of state, the critical temperature 
Tc was calculated to vary in a range of around 15 to 20 
MeV for nuclear matter HE El El- (One of the latest 
investigations, based on an improved Fisher's model J23| 
indicated Tc — 6.7 ± 0.2MeV for finite nuclear systems. 
This value is source of controversy, e.g. (lH^ll)- Dur- 
ing the reaction process, the system explores different 
regions of the phase diagram. Since at relativistic en- 



ergies the collision is related to short wavelengths, the 
hot remnant should reach high positive values of pressure 
P due to thermal energy (rather than mechanical com- 
pression, characteristic of Fermi-energy collisions) with- 
out deviating sensibly from the initial density po- It is 
commonly assumed that at this stage the system is still 
not thermalised and it undergoes expansion in order to 
attain equilibrium (There exist also opposite interpreta- 
tions assuming thermalization already before expansion 
and a successive "Big-Bang-Hke" expansion out of equi- 
librium (^). If the initial pressure is high enough, the 
subsequent expansion could lead to rather low densities, 
and the system, after dissipating the incoming momen- 
tum, could reach a point belonging to the spinodal region. 
Due to the inverse relation between pressure and density 
dP/dp < 0, this region is unstable, and density fluctu- 
ations are magnifled. The nucleus breaks apart due to 
spinodal instability. The system disassembles also due 
to Coulomb instability. The inclusion of the long-range 
Coulomb interaction in the equation of state was intro- 
duced by Levit and Bonche [23|, with the result that 
the solution of the coexistence equation vanishes above a 
"limiting temperature" Tum, in general much lower than 
Tc, depending on the conditions taken for the calcula- 
tion (see also [sl, l^^l)- Density fluctuations reflect a 
continuous evolution of the size and number of nuclear 
droplets from a conflguration to another yJI- If the av- 
erage mutual distance among the nucleons exceeds the 
strong nuclear interaction range (i.e. about \J < (t„ > /tt, 
where < cr„ > is the average nucleon-nucleon collision 
cross section), the break-up conflguration "freezes" and 
the formed nuclei and nucleons fly away freely, all carry- 
ing signatures of the so-called freeze-out temperature of 
their common source. From comparing results from dif- 
ferent experimental approaches e.g. (iM |2J, l3flL l8q this 
temperature is found to be restricted to a range of 5 to 6 
MeV (corresponding to a range of excitation energy per 
nucleon around 2.5 to 3.5 MeV), quite independently of 
the reaction. This flnding, not directly compatible with 
the phase diagram of ideal nuclear matter even suggested 
to search for a "characteristic temperature" of fragmen- 
tation [23. The break-up conflguration at freeze-out is 
expected to reflect the excitation energy of the system. 
The dense phase of highly heated systems should have 
the aspect of an ensemble of copious almost-equal-size 
light fragments. At reduced excitation, just sufHcient 
for attaining the freeze-out, the break-up partition might 
evolve to more asymmetric conflgurations, where the for- 
mation of a heavy fragment close to the mass of the hot 
remnant is accompanied by one or more light fragments 
and clusters. As an extreme, this conflguration might 
even reduce to a binary asymmetric decay. In the case 
of a very asymmetric split of the system, the partition 
multiplicity has minor influence on the kinematics of the 
light ejectiles. The emission of light particles populates 
spherical shells in velocity space and can not be easily dis- 
tinguished by the kinematics from a binary decay when 
large mass-asymmetries characterize the partition. A bi- 



nary or binary-like decay issued from a break-up con- 
flguration is a "fast" process. Compared to asymmetric 
flssion, asymmetric break-up decays should result in a 
similar U-shape of the mass spectra of the yields. On 
the other hand, break-up decays should be reflected in 
the higher magnitude of the yields, and in the emission 
kinematics that, still mostly governed by the Coulomb 
fleld, should exhibit an additional contribution due to 
the eventual expansion of the source. 



B. Measurement of light-fragment properties 

Great part of the information on light-particle emission 
at high excitation energies was collected in 4-7r-type ex- 
periments, suited for measuring the multiplicity and the 
correlations of intermediate-mass fragments [3^133. l4(l|. 
Still, the measurement of correlations and the linear- 
momentum-transfer was the basis for pursuing intense 
researches on the transition from the formation of com- 
pound nuclei to multifragmentation jij |42| • 

In this work, we discuss additional results derived from 
new inclusive measurements of the reactions ^^Fe+p and 
^^Fe+"^*Ti at 1 ^ GeV, effectuated in inverse kinematics 
with the FRagment Separator (FRS) at GSI (Darm- 
stadt). The experimental set-up was not intended to 
measure multiplicity and correlations, but to provide for- 
mation cross sections and high-resolution velocity spectra 
for isotopically identifled projectile-like residues. The ex- 
citation of the ^^Fe-l-p system consists of purely thermal 
energy, and it is just high enough to approach the condi- 
tions for the onset of multifragmentation. On the basis 
of these data we search for the properties of the early ap- 
pearance of break-up events and their competition with 
compound-nucleus emission. The system ^^Fe+^'^^Ti is 
compatible with an abrasion picture. The excitation en- 
ergy deposited in the projectile spectator, still mostly of 
thermal nature, establishes the dominance of multifrag- 
mentation in the decay process. We will especially dis- 
cuss the differences in the kinematics of light-fragment 
emission in the two systems, conditioned by two differ- 
ent levels of excitation magnitude. 



II. EXPERIMENT AND ANALYSIS 
PROCEDURE 

The experiment was performed at GSI (Darmstadt). 
A primary beam of ^^Fe was delivered by the heavy-ion 
synchrotron SIS at an energy of 1 A GeV. The target 
was constituted of liquid hydrogen (with a thickness of 
87.3 mg/cm^) contained in a cryostat with thin titanium 
windows (36.3 mg/cm^ in total), wrapped in thin Mylar 
foils (C5H4O2, total thickness: 8.3 mg/cm^) for thermal 
insulation. In the target area, other layers of matter in- 
tersected the ion-beam: the accelerator-vacuum window 
of titanium (4.5 mg/cm^) and the beam-current moni- 
tor composed of aluminium foils (8.9 mg/cm^). In or- 



der to disentangle the production and the physical re- 
sults related to the interaction with hydrogen from the 
contribution associated to the other materials, the whole 
experimental runs were repeated in identical conditions, 
after replacing the target by titanium foils having the 
same thickness of the cryostat windows and wrapped in 
Mylar foils having the same thickness of the cryostat in- 
sulation. This procedure did not only determine the dis- 
turbing contributions in the measurement of the ^^Fe-|~p 
system, but it also provided additional experimental data 
on an other reaction system. With some arbitrariness 
we name "titanium target (natfjiy' ^j^g ensemble of the 
titanium foils replacing the cryostat window, the My- 
lar wrapping, the accelerator-vacuum window and the 
beam-current monitor. Unfortunately, the measurement 
of the ^^Fe+"^*Ti system accounts also for non-titanium 
nuclei, the pollution of which corresponds to their por- 
tion in the total number of target nuclei per area and is 
equal to 25.9% (Al) + 7.2% (Mylar) = 33.1%. It should 
be remarked that these components are not placed at 
the same distance from the entrance of the spectrometer. 
Fragments produced in the beam-current monitor or in 
the accelerator-vacuum window could have lower proba- 
bilities to be registered in the experiment since the angu- 
lar acceptance is reduced by factors of 0.33 and 0.25, re- 
spectively, compared to products from the titanium foils 
replacing the cryostat. Henceforth, we refer to the liquid 
hydrogen as "proton target (p)". In this case no polluting 
contributions are included in the final results. 



experiment are shown. Events are ordered according to 
the measured Z and AjZ so as to obtain an isotopic 
identification plot. 



B. Longitudinal velocities 

The measurement of the time of flight is precise enough 
for an accurate identification of the mass of the frag- 
ments. Nevertheless, mainly due to the resolution and 
additionally due to a slight dependence on the trajec- 
tory J4J| it is not suited for a fine measurement of the 
velocities of the fragments. On the other hand, once an 
isotope is identified in mass and charge, a much more 
precise measurement of the velocity is obtained directly 
from the magnetic rigidity of the particle 
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In this case, the precision of /?7 depends only on Bp, 
that has a relative uncertainty of 5 ■ 10~^ (FWHM) for 
individual reaction products. The absolute caHbration of 
the defiection in the magnet in terms of magnetic rigidity 
Bp is performed at the beginning of the experiment with 
a dedicated calibration run using the primary beam. 

Since one single magnetic configuration of the FRS se- 
lects only a Bp range of about ± 1.5%, several overlap- 
ping runs have been repeated imposing different mag- 
netic fields. While for the heavy residues close to the 



A. Nuclide identification 

The reaction products were analysed with the high- 
resolution magnetic spectrometer FRS, constituted of 
four dipoles set in an achromatic mode. The horizon- 
tal positions of the reaction products were deduced from 
two plastic scintillators, the first positioned in a disper- 
sive plane after the first two dipoles, the second installed 
in the final achromatic plane. From the time of flight 
(TOF), evaluated between the two scintillators, the rel- 
ative velocity (3 and the 7-factor were obtained. The 
particle charge q was measured with an ionisation cham- 
ber placed in front of the achromatic plane. Since the 
reaction products were fully stripped, the nuclear charge 
Z = q was deduced directly. The mass A was deduced 
from the time of flight and the magnetic rigidity of the 
particles according to the relation 
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where Bp is the magnetic rigidity of a particle, c the 
velocity of light, e the elementary charge, mo the nuclear 
mass unit, Sm ~ dM/A the mass excess per nucleon. For 
the purpose of the isotopic identiflcation, the variation of 
Sm with A/Z can be neglected, and a linear variation of 
A/Z as a function of Bp/Pj can be assumed. In fig. |2l 
the raw data collecting all the events measured in the 



Z = 20 




FIG. 2: Experimental isotopic resolution. The isotopes are 
grouped in chains. One chain collects nuclei having the same 
value of the difference N — Z. The straight chain corresponds 
to A/Z = 2 or N - Z = 0. On the right of this chain the 
isotopes are neutron rich, on the left they are proton rich. 



projectile one or few settings were sufficient to cover the 
whole velocity spectrum, the light fragments often re- 
quired more than ten runs. The Bp scanning of ^Li, 
produced in the interaction with the target of liquid hy- 
drogen enclosed in the cryostat constitutes the diagram 
(a) of fig. 13 each segment of the spectrum is obtained 
from a different scaling of the set of magnetic fields of 
the FRS. In order to obtain consistent weightings, the 
counts of the different measurements were normalized to 
the same beam dose. For each magnetic scaling, this 
normalization was obtained by dividing the correspond- 
ing segment of the spectrum by the number of projectiles 
that hit the target during the corresponding run. The im- 
pinging projectiles were counted with the beam-current 
monitor. The renormalized yields are shown in the dia- 
gram (b) of fig. 13 We should note that the spectrometer 
accepts only the fragments emitted in a cone of about 15 
mr around the beam-axis in the laboratory frame, when 
the reaction occurs in the hydrogen-target position. As 
a consequence, a light residue like, for example, ^Li, gen- 
erated in a collision at a beam energy of 1 ^ GeV can be 
detected only if emitted with small transverse momen- 
tum. The experimental spectrum represents the part of 
the density distribution in the velocity space selected by 
the angular acceptance of the spectrometer, projected on 
the longitudinal axis. Unfortunately, the angular accep- 
tance depends on the magnetic rigidity of the particles. 
As pointed out in the work J4a| , for a given set-up of the 
spectrometer, the more the intersection of the trajectory 
of a particle with the dispersive or the achromatic planes 
is displaced from the centers, the lower is the acceptance 
angle of the FRS. The effect appears in the curved sides 
of each single segment, with the result of disturbing the 
overall structure of the Bp scanning. This distortion, 
seen in the spectrum of the plot (b) of fig. |3 can be suc- 
cessfully corrected by means of ion-optical calculations 
that fix the dependence of the angular transmission on 
the trajectory. The calculation of the ratio of the trans- 
mission T relative to its maximum value is presented in 
fig. E The corrected spectrum, seen in the plot (c) of 
fig. 13 is the result of scaling up the yields of the spec- 
trum by the factor T-^ax/T. We also changed from a 
/?7 spectrum to a longitudinal-velocity spectrum and, to 
simplify the analysis, the reference frame was changed 
from the laboratory to the beam frame. On the aver- 
age, the projectile interacts in the middle of the target. 
Therefore, we take into account the slowing down of 
^^Fe in the first half of the target, as represented in the 
upper diagram of fig. [3 We also consider that the frag- 
ments slowed down in the remaining half of the target 
and, therefore, were emitted at higher velocity than the 
one we observed. The analysis so far illustrated was re- 
peated for all the isotopes produced in the interaction 
with the target of liquid hydrogen enclosed in the cryo- 
stat. Successively, the same procedure was applied to 
the corresponding isotopes produced in the interaction 
with the "'^^Ti target. As all spectra are normalized to 
the same beam dose, by subtracting the velocity spectra 
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FIG. 3: Four steps of the analysis procedure to obtain the 
observed velocity spectrum of "^Li emitted in the reaction 
^^Fe+p. (a) Raw spectra of counts as a function of f3j in 
the laboratory frame. Each segment results from a different 
scaling of the magnetic fields of the FRS. One segment asso- 
ciated to the same magnetic scaling is marked with hatched 
areas in this plot and in the two following ones. Arrows de- 
limit the scanned fiy range, (b) Yields normalized to the same 
beam dose, (c) Elimination of the angular-transmission dis- 
tortion. Spectrum as a function of the longitudinal velocity 
in the beam frame wji . The broad Gaussian-like hatched area 
indicates the contributions from non-hydrogen nuclei, (d) All 
components of the spectrum are composed together averaging 
overlapping points. Contributions from non-hydrogen-nuclei 
were suppressed. The spectrum was divided by the number 
of nuclei per area of the liquid- hydrogen target. Statistical 
uncertainties and a fit to the data are shown. 



of the residues produced in 56pg^natrp| (indicated by the 
hatched area in the plot (c) of fig. O from those of the 
corresponding isotopes produced in the target of hydro- 
gen stored in the cryostat, we could obtain the measured 
velocity distributions for the reaction with the liquid hy- 
drogen. The resulting yields are unambiguously disentan- 
gled from any disturbing contributions produced by other 
material present in the target area. Finally, the velocity 
spectra obtained for the ^^Fe-|-p system were divided by 
the number of nuclei per area of the proton target. The 
resulting spectrum is shown in the diagram (d) of fig. |3 



In the case of the ^^Fe+"'**Ti system, we should consider 
that the target is constituted of three components, the 
titanium foils replacing the cryostat, the beam-current 
monitor and the accelerator-vacuum window, having a 
number of nuclei per area equal to ng, ni, and 712, re- 
spectively. We should also recall that these components 
are placed at different distances from the entrance of the 
spectrometer and are subjected to different values of the 
angular acceptance, that is about a^ ~ 15.8 mr, ai = 9 
mr, and 0:2 = 7.8 mr, for the layers uq, ni, and n2, re- 
spectively. Thus, the cross sections given in this work for 
the "titanium" target are calculated using a target com- 
position where the different layers are weighted by the 
corresponding estimated transmission values T, assum- 
ing identical production cross sections in the different 
target components. In particular, the velocity spectra 
obtained for the ^^Fe+"^*Ti system should be divided by 
the quantity noT{ao) + niT{ai) + n2T{a2). 

The experimental data are already complete enough to 
let us recognise an important signature of the Coulomb 
repulsion: the double-humped spectrum reveals that the 
velocity of ^Li nuclei emitted at small angles has two 
components: one appreciably higher and one appreciably 
lower than the beam. According to the references J4a, 
|43,|4^, where similar structures have been observed for 
fission fragments, we may connect the double-humped 
spectrum to the action of the Coulomb field of a heavy 
partner in the emission process. 

Once changed to longitudinal velocities in the beam- 
frame Wm, the shift of the barycenter of the spectrum 
with respect to zero is equal to the mean reaction recoil 
< Um >. Also this quantity, studied in the lower dia- 
gram of fig. carries a valuable information about the 
reaction mechanism, and it can be related to the friction 
suffered by the projectile in the collision, according to a 
given impact parameter [43|. Due to the Hmited angu- 
lar acceptance of the FRS which favours the detection of 
heavier nuclei, a depletion of the statistics for the mea- 
surement of the lightest nuclei is expected when, as in 
this case, the light fragments are measured together with 
the heaviest in the same magnetic setting. Due to these 
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FIG. 4: Transmission of the FRS as a function of the posi- 
tions in the dispersive and achromatic planes, relative to its 
maximum value. Numerical values are taken from ref. lifill. 
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FIG. 5: Top. Definition of the beam frame and of the center- 
of-mass frame of the emitting source with respect to the labo- 
ratory frame. The diagram corresponds to realistic conditions 
of the present experiment for ^Li. The solid lines describe the 
slowing down of the beam and of the centroid of the veloc- 
ity spectrum of ^Li in traversing the target. Bottom. Mean 
longitudinal recoil velocities in the beam frame < uS > of 
the reaction residues compared with the systematics of Mor- 
rissey |4g| (solid line); only isotopes with sufficient statistics 
and entirely measured velocity spectra are considered. 



problems, the mean velocities of the light residues can 
only be determined with relatively large uncertainties. 
With these large uncertainties, the information from the 
mean velocities could not be exploited. Although these 
mean velocities also enter into the evaluation of the cross 
sections, the uncertainties they introduce are compara- 
ble to those from other sources. We preferred to deduce 
the mean recoil velocities of lithium, beryllium, boron 
and carbon by extrapolation from the systematics of the 
data relative to the ensemble of the heavier residues. 



III. 



RESULTS 



A. Nuclide cross sections 

When a fragment is emitted with a large absolute ve- 
locity f = |f/| in the center of mass, not all the angles 
of the corresponding velocity vector v are selected by 
the finite angular acceptance of the spectrometer. As a 
result of the data analysis detailed in the previous sec- 
tion, we obtain the measurement of the apparent cross 
section dr(w|[)/du|| as a function of the longitudinal ve- 
locity W|| . This observed cross section differs from the real 
cross section due to the angular acceptance. The detec- 



tion of a particle depends on the perpendicular velocity 
v± = y/v'^ — Wn in the center-of-mass frame, the angle of 

rotation around the beam direction (p, and the velocity u 
of the center of mass with respect to the laboratory. The 
dependence on ip comes about because the beam pipe in- 
side the quadrupoles is not cyHndrical. To reconstruct 
the full velocity distribution, independent of the angu- 
lar acceptance of the spectrometer, an assumption on 
the angular distribution is necessary. It was concluded 
from experiments, to which the full angular range was 
accessible, that the data are in satisfactory agreement 
with an isotropic emission (see, for example, the treat- 
ment of "Moving source analysis" presented in j^). This 
assumption has been corroborated by a vast collection 
of data for reactions of very different nature. Isotropic 
emission has been observed either for lowly excited fis- 
sioning systems 0, or even for very highly excited nu- 
clei under goin g expansive flow in thermal multifragmen- 
tation (lftll7l|. At least the ^^Fe-|-p system can be safely 
included in this range. Slightly less justified is the as- 
sumption for ^^Fe+^'^'Ti, since some effects of dynamical 
multifragmentation could disturb the isotropy. 

Thus, if we assume isotropic particle emission in the 
center-of-mass frame, the variation of the cross section 
cr{v), as a function of the absolute velocity v, is related 
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FIG. 6: Experimental isotopic production cross sections of 
some light elements for the reactions *^Fe+p and ^^Fe+'^'^'Ti 
at 1 A GeV. The cross sections related to the latter system 
are scaled of a factor 0.07. Numerical values are collected in 
table HI 




FIG. 7: Isotopic production cross sections shown on a chart 
of the nuclides for the reactions '^'^Fe+p and '^'^Fe+"''*Ti at 
1 A GeV. The values for ^"Mn in ^'^Fe+p and for ^^Cr in 
^^Fe+^'^^Ti were obtained from systematics. 



to the variation of the apparent cross section X(w|j) as a 
function of vi\ by the equation (see appendix I A l|l : 
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where v± = v±{(p,u) is the highest value of v± selected 
by the angular acceptance of the spectrometer, and ip is 
the rotation angle around the beam direction. As demon- 
strated in the appendix lA 21 it is possible to reverse the 
relation (pj and extract da{v)/dv by the straightforward 
solution of a system of geometric relations. The forma- 
tion cross sections are directly obtained by integration of 
da{v)/dv. In the appendix^l tableUcollects the isotopic 
cross sections for the production of light residues, from 
lithium up to oxygen, measured in this work for the re- 
action ^""Fe-l-p and ^'^Fe+"''*Ti. The distributions of the 
formation cross sections evaluated for the two systems 
^^Fe+p and sepg^natrj,- ^^u^ Q^y ^^.g presented in fig.EI 
for different elements as a function of the neutron number 
and in fig.d on the chart of the nuclides. The extension 
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FIG. 8: Experimental production cross sections as a function 
of the mass number. The statistical uncertainties are lower 
than 10%. The systematic uncertainties evolve from 10% for 
the heavy residues close to the projectile to 20% for the light 
fragments. 



in the two systems. On the contrary, the particularity of 
the proton-induced spallation compared to the titanium- 
induced fragmentation arises strikingly when the kine- 
matics of the light-particle emission is investigated. The 
density of velocity vectors z/ in a plane containing the 
beam axis is presented in fig. El (As observed in the fol- 
lowing, this presentation is equivalent to invariant-cross- 
section plots.) Thus, when we compare the ^^Fe-|-p sys- 
tem to the ^^Fe+"^*Ti system on the basis of the recoil 
velocity, we find that substantially distinct mechanisms 
should be involved in the light-fragment emission in the 
two systems. In the fragmentation reaction induced by 
the titanium target, all the residues are emitted accord- 
ing to a bell-shape velocity spectrum. A long sequential 
decay would produce this kind of shape; in this process, 
neutrons, protons and clusters are in fact emitted with 
different angles with equal probability. Nevertheless, 
due to the high excitation of the hot fragments gained 
in such a violent collision, and due the exponential in- 
crease of the cross section of the light residues with the 
mass loss, we are in favour of a multifragmentation pic- 
ture to depict the dominant deexcitation process. In this 



of the production appears rather similar and, in partic- 
ular, despite the expected difference in excitation energy 
reachable in the coUisions with the two different targets, 
the cross-section distributions of the residues of the two 
reactions do not manifest drastic differences in their fea- 
tures. A more quantitative revelation of this similarity is 
presented in fig.|Hl where the mass distributions are com- 
pared. The difference in the shape of the mass spectra 
is significant only for the intermediate masses: the cross 
section of the residues of ^^Fe+p decreases from A = 30 
to A = 18 by about one order of magnitude, while we 
observe only a slight decrease by about a factor of two 
for ^^Fe+"^*Ti. The data reveal that higher excitation 
energy introduced by the interaction with titanium, with 
respect to proton-induced spallation, results in decreas- 
ing the slope of the mass-spectra in the IMF-range and 
depleting the cross section for heavy residues in favour 
of an enhanced production of light fragments. However, 
the portion of the mass spectra corresponding to light- 
particle emission follows a very similar exponential slope 
for both systems. 
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B. 



Velocities 



In the previous section, the full velocity distributions 
were reconstructed from the data and employed to obtain 
the residue formation cross sections for the reactions in- 
duced by the proton and titanium target, respectively. 
Though, the cross sections did not yield any unambigu- 
ous distinction between the two reactions that, indeed, 
should result into rather different deexcitation pictures 
on the basis of the different thermal excitations reached 
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beam frame {v,. 



representing the distribution on a plane 



containing the beam axis. The solid lines denote the angular 
acceptance of the spectrometer. 
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FIG. 10: Velocity spectra of light residues produced in Fe+"''*Ti (upper diagram), and in ^ Fe+p at 1 y4 GeV (lower diagram), 
ordered on a nuclear chart. They are represented as a function of the velocity in the beam direction in the beam frame Wy . 
Crosses and points indicate measured spectra dX ^ /Avu and dl^ jAvu, respectively, defined according to eq. (Ol, and normalized 
to the unit. They represent all fragments transmitted through the FRS. Reconstructed velocity spectra crj' and crj?, defined 
according to eq. Q and normalized to the unit are marked with dashed and solid lines, respectively. In the lower diagram, the 
reconstructed spectra for Li, B and C emitted from Fe+"^'Ti are superimposed as dashed lines for comparison. 



10 



case, the hot source is expected to undergo a fast ex- 
pansion and successively form several fragments. In this 
scenario [ij| , the emission velocity of a light residue could 
vary largely according to different parameters: the par- 
titioning in the multifragmentation event, the expansion 
of the source before the break-up phase, and the position 
where the tracked fragment is formed with respect to the 
other fragments. Also this process would result in a ve- 
locity spectrum with a bell shape centerd at the mean 
recoil velocity, equal to the shape we observe. 

On the contrary, when light fragments originate from 
the ^^Fe+p system, the reaction dynamics leads to the 
population of one most probable emission shell in the ve- 
locity space, around the center of mass. This is the case 
of ^Li, as shown in fig. El Only a forward and a back- 
ward portion of the emission shell could be measured as 
selected by the conical cut that the spectrometer deter- 
mined: this fully explains the double-humped velocity 
spectra shown in fig. El 

The velocity distributions of the light fragments gener- 
ated in the proton-induced reaction carry the unambigu- 
ous signature of a strong Coulomb repulsion in the emis- 
sion process. This observation evidently excludes that 
the light fragments could be the final residues of a long 
sequential evaporation chain. The strong Coulomb com- 
ponent in the emission process rather reflects the domi- 
nating influence of a very asymmetric split of the source. 

We can now reduce the representation of the recoil- 
velocity distribution a{v) to one dimension, selecting 
only those velocities v aligned in the beam direction, and 
occupying only abscissae in the plots of flg. El 

Due to our assumption of isotropy, we can deflne ra- 
dial velocity distributions dividing the differential cross 
section dcr(w)/du associated to a given velocity v in the 
center of mass by the spherical surface of radius v. 



dV 



1 da 

Anv^ dv 



(4) 



It should be remarked that either in the reference of the 
center of mass or in the projectile frame, 7 is close to 
the unit and consequently Cr (w) is directly related to the 
invariant cross section <ji{v). Indicating m = jm the 
mass of the particle, p = 77? its momentum and E its total 
energy in the center of mass frame (or in the projectile 
frame), we obtain the equality: 
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Also the planar cuts in velocity space {v^^,v']^) of flg. El 
are equivalent to invariant-cross-section plots |^ . 

As a technical remark, the advantage of inverse 
kinematics compared to direct-kinematics experiments 
should be pointed out. The registration of emission ve- 
locities close to the velocity of the center of mass of the 
hot remnant are not prevented by any energy threshold. 
Thus, only in inverse kinematics we can clearly appreciate 
the gradual transition from a chaotic-dominated process. 



reflected in Gaussian-like invariant-cross-section spectra, 
to a Coulomb- (or eventually expansion-) dominated pro- 
cess, producing a hollow around the center of mass. This 
characteristic signature we exploit resembles the investi- 
gation of relative velocity correlations between two frag- 
ments j^ in full-acceptance experiments for analysing 
decay times. In that case, the probability to detect two 
almost simultaneously emitted particles in space with 
small differences in direction is suppressed due to the 
mutual Coulomb interaction. 

A systematic study of the spectra of lithium, beryl- 
lium, boron and carbon is presented in flg. ^1 and com- 
pared with the observed velocity distributions. In the 
^^Fe+"^*Ti reaction, all spectra show a bell shape. In 
the ^^Fe+p spallation, the double-humped distribution 
appears clearly for isotopes with mass lower than twelve 
units. The shape of the velocity spectra depends mostly 
on the mass rather than on the charge, and chains of 
isotopes belonging to the same elements show a transi- 
tion from a bell shape toward a double-humped spectrum 
with decreasing mass. This transition is not always grad- 
ual but, as revealed by the neighbouring ^^C and ^^C in 
the lower panel of flg. El sometimes seems to be rather 
abrupt. 



IV. 



DISCUSSION 



A. Systematics of kinematical features 

On the basis of the ensemble of experimental data on 
the production cross section and on the emission veloc- 
ity of the residues, we devote this section to discuss the 
reaction phenomenology. In accordance to the vast lit- 
erature dedicated to ion-ion fragmentation (explored in 
the reviews [Hiil), we can safely relate the ^'^Fe+^'^'Ti 
reaction to the formation of highly excited systems, the 
decay of which is commonly interpreted as a multi-body 
instantaneous disassembly. In the following, we will refer 
to the ^^Fe+"^'Ti collision as a guideline for comparing 
to a fragmentation scenario. We will rather concentrate 
on the reaction mechanism of the proton-induced colli- 
sion, which points to competitive types of decay at the 
same time. 

From a flrst analysis, the main kinematical characterit- 
ics of the ^^Fe-|-p system, recalling a strong Coulomb re- 
pulsion, evocate a binary decay process. 

On the other hand, other relevant features, like the 
high production yields for both light and about half- 
projectile-mass residues, could evocate the character of 
a fast decay, in line with the scenario of a sudden disas- 
sembly of the source depicted for the ^^Fe+^'^^Ti system. 

Our flrst attempt will be to test the pertinence of the 
experimental data on the emission velocities with a gen- 
eral systematics. Afterwards, we will discuss the intri- 
cacy of the several possible contributions to the spec- 
tral shape of the kinetic-energy distributions, and the 
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difficulty to extract insight on the excitation energies in- 
volved in the reaction directly from the measurement. 



1. Absolute-velocity spectra 

A recurrent analysis of the Coulomb-repulsion aspects 
is the comparison of the distribution of absolute veloc- 
ities of outgoing fragments v = \v\ (where v is the cor- 
responding velocity vector in the center of mass of the 
hot remnant) with the systematics of total kinetic en- 
ergy released in fission. We intend to follow this ap- 
proach (e.g. j^) to test the compatibility of the light- 
fragment emission in ^^Fe-l-p with an asymmetric-fission 
picture. The FRagment Separator is particularly effi- 
cient in measuring recoil velocities, because the magnetic 
rigidity of the residues is known with high precision (see 
section ITTb)| . Indeed, the identity of the mother- nuclei 
is hidden in the complexity of the interaction processes 
related to high-energy collisions, like the intra-nuclear 
cascade and some evaporation events prior to the binary 
decay. The present new data are especially significant as 
they are the first measurement of the velocities of frag- 
ments issued of proton-induced splits of iron-like nuclei. 
On the other hand, fission velocities of residues of fight 
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FIG. 11: Left panel Mean absolute velocities in the reference 
frame of the center of mass of the hot remnant, measured 
for residues of the ^^Fe-fp system (open circles) and deduced 
from the systematics of Tavares and Terranova [^ (hatched 
bands). The width of the hatched areas results from the range 
of the possible mother nucleus from *^Ti (lower values) to 
^ Fe (higher values). In the insert, data points on the total 
kinetic energy released in a symmetric split of nuclei close 
to iron, measured by Grotowski et al. [54l | are compared to 
the systematics of Viola [5fil | (dot-dashed line) and with the 
systematics of Tavares and Terranova [^ (solid line). Right 
panel Measured absolute-velocity spectra for the residues Li, 
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^°B, and '^^C produced in the '^"Fe-fp system. The arrows 
indicate the values obtained by the systematics of Tavares 
and Terranova. 



nuclei have been widely investigated in fusion-fission ex- 
periments |j|, with the advantage of excluding most of 
the ambiguities on the identification of the fissioning nu- 
cleus. Data on symmetric fission of nuclei close to iron, 
formed in fusion reactions were published by Grotowski 
et al. 01 and where the basis for the revised kinetic- 
energy- release systematics of Viola [^. This systemat- 
ics establishes a linear dependence of the most probable 
total kinetic energy El. released in a symmetric fission to 
the quantity Z'^/A^^'^, evaluated for the mother nucleus: 



El 



aZ^A 



1/3 



(6) 



where A and Z identify the fissioning nucleus and a 
and b are parameters fitted to the experimental data 
(a = (0.1189 ± 0.0011) MeV, b = (7.3 ± 1.5) MeV). 
More recently, new data obtained for the binary split 
of even lighter nuclei than iron inspired Tavares and Ter- 
ranova j^ to revisit the systematics of Viola once more. 
The new systematics is close to the systematics of Viola 
for heavy nuclei down to Z'^/A^'^ w 200. As shown in 
the insert of fig. E] iron-like nuclei constitute a turn- 
ing point: for lower masses the function changes slope, 
so that the total kinetic energy released vanishes for Z 
approaching 0. As anticipated by Viola [^, the expec- 
tation for a slope change around iron results by the effect 
of diffuseness of light nuclei in disturbing the formation 
of the neck, in the liquid-drop picture. The following 
relation was deduced: 
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where 



b and c are fitting parameters (a = 



9.39MeV-\ 6 = -58.6MeV-i, c = 226MeV-i). 

Since the systematics is valid for symmetric splits only, 
a term should be added to extrapolate to asymmet- 
ric splits, when two fragments are formed with masses 
mi, 1712, mass numbers Ai,A2, and charges Z\,Z2, re- 
spectively. Following the hypothesis of non-deformed 
spheres at contact (as also imposed in jSol), the Coulomb 
potential is proportional to the product of the charges of 
the fission fragments Z\Z2, divided by the distance of 

1/3 1/3 

their centers, which varies with A^ ^ A^ . The con- 
version from the symmetric to asymmetric configuration 
is therefore : 



(8) 



It should be remarked that the possible presence of a neck 
is not included in this simple relation that, therefore, is 
a good approximation for light systems only. From the 
momentum conservation and the introduction of the re- 
duced mass /z = miTO2/(TOi -I- m^), we can relate the 
total kinetic energy to the velocity v\ of the fragment A\ 
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by the relation E^ 



form of El in the relations 



m\v1/2^. Introducing the latter 
l|Hl, and substituting the to- 
tal kinetic energy released in symmetric fission with the 
corresponding value given by the systematics E^^^^, we 
obtain the conversion 
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Following the strategy of previous publications, e.g. [^, 
for the same light residues we compare the centroids of 
the measured absolute velocity spectra to the predicted 
velocities in fission events; the latter are deduced from 
the systematics of total kinetic energy released in fission 
by applying the relation Jsj. In the right side of fig.lTTl 
we observe that, while the most probable absolute ve- 
locity does not diverge considerably from the systemat- 
ics (assumed for ^^Fe as mother nuclei), the spectra of 
lighter fragments exhibit a long exponential tail for very 
high velocities. As a consequence of the asymmetry of 
the absolute velocity spectra with respect to a Gaussian 
distribution, the experimental centroids lie above the fis- 
sion systematics, as shown in the left side of fig. El The 
hatched bands represent the range in velocity due to an 
assumed variation of the mother nucleus from "^^Ti (lower 
velocities) up to the projectile (higher velocities). In pre- 
vious works (e.g. |l6, 53, aj|) such tails to very high ve- 
locities, reflected in the divergence from the systematics, 
were related to the emission from an expanding system 
in its initial expansion stage. 



2. Kinetic- energy spectra 

Directly obtained from the absolute-velocity spectra, 
the distributions of kinetic energy Ek offer another rep- 
resentation of the kinematics, where some more classic 
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FIG. 12: Kinetic-energy spectra in the center of mass of the 
emitting source obtained from the reconstructed experimental 
velocity distributions in the case of emission of "^Li and ^^C, 
respectively. The spectra are compared for the ^^Fe-|-p system 
(solid lines) and for the ^ Fe+"''*Ti system (dashed lines). 
All spectra are normalized to the same area. The smooth 
distributions result from a spline fit procedure to the data. 



features could be searched for. In flg. 1121 similarities and 
differences in kinetic-energy spectra associated with pro- 
ton and titanium target nuclei are illustrated. The tails 
to high emission velocities (flg. E] right panel) lead to 
long tails in the kinetic-energy spectra and characterize 
both systems. We interpret it as a general indication that 
the collision generated very high excitation energy in the 
system. It would be tempting to even deduce the ther- 
mal properties of the system. In this case, with particular 
concern for the ^^Fe+p system, we could draw assump- 
tions on the probability for break-up channels. Unfortu- 
nately, even if in some studies the nuclear tem per ature 
was deduced from the inverse slope parameter [23|, the 
mixing up of several effects in the observed kinematics 
yields serious ambiguities in the extraction of thermal 
properties of the source. We can list at least eight of the 
combined effects describing the observed spectral shapes. 

1. The presence of a Coulomb barrier results in the 
deviation of the spectral shape from a Maxwell- 
Boltzman distribution (the maximum moves to 
higher values). 

2. The transmission through the barrier is ruled by 
a Fermi function with an inflection point at the 
barrier and not by a discontinuous step function. 
This effect introduces a widening of the spectrum. 

3. As a result of the initial stage of the collision, an 
ensemble of several possible sources with different 
Z and A are related to different Coulomb barri- 
ers. The folding of different Coulomb barrier peaks 
results in a broader hump. 

4. If emitted nuclei undergo further evaporation 
events, the spectrum widens. 

5. The temperature of the hot source acts on the 
recoil momenta of the emitted fragments. If at 
least major disturbing effects like the variation of 
the emitting source, the Fermi momentum in the 
hot fragmenting nucleus, described below, and the 
transmission through the barrier were negligible, it 
would be possible to deduce the temperature of the 
equilibrated fragmenting system from the inverse 
slope parameter fltted to the tail of the high side 
of the energy spectrum of the residues. 

6. The Fermi momentum of particles removed in the 
collision with protons or abraded in the interac- 
tion with the titanium target produces a momen- 
tum spread that could be evaluated according to 
Goldhaber's formaHsm [531. 



PF 



A-1 
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where A is the mass of the hot remnant, Ai is the 
mass of the emitted cluster and ap is the Fermi- 
momentum spread. The momentum spread deriv- 
ing from the Fermi-momentum spread produces a 
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distribution of momenta of the center of mass of the 
remnants in the projectile frame. In deducing the 
energy spectra of the residues in the frame of the 
center of mass of the remnant, the spread related 
to the Fermi-momentum could not be eliminated as 
the mass of the remnants are unknown. As a result, 
the Fermi-momentum contributes both to widening 
the spectrum and incrementing the tail for high en- 
ergies. Quoting from Goldhaber [^, when a ther- 
maHsed system with a temperature T and mass A 
emits a cluster of mass Ai , the momentum spread 
of the fragment spectrum is 
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where mo is the nuclear mass unit and k is Boltz- 
mann's constant. The momentum spread ap,, re- 
lated to the Fermi momentum adds to the momen- 
tum spread induced by the reaction. This means 
that, just reversing the previous relation, the ad- 
ditional contribution to the temperature related to 
the Fermi momentum is equal to the apparent tem- 
perature 
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As it was remarked in early studies [^, the extrac- 
tion of the nuclear temperature from the measured 
energy spectra of the residues is therefore a danger- 
ous procedure (a recent discussion of the problem 
of the Fermi motion is presented in jSfll)- 

7. Multifragmentation events could be accompanied 
by the expansion of the nuclear system. Nuclei 
emitted in the initial instant of the expansion would 
populate the high-energy tail of the s pec trum. This 



is the case for very excited systems 



i pec t 



8. The multiplicity of intermediate-mass fragments 
simultaneously emitted might be reflected in the 
maxima. According to previous investigations |62 |. 
a drop in the maximum energy of the outgo- 
ing fragments in a simultaneous disintegration of 
the source indicates higher average multiplicity of 
intermediate-mass fragments: this is related to the 
larger number of participants in the redistribution 
of the kinetic energy. 

The last of the enumerated contributions to the energy 
spectra is evident in fig. El In the proton-induced 
collision, the position of the maximum corresponds to 
larger kinetic energy than in the case of the titanium 
target. This might be related to higher multipHcity of 
intermediate-mass fragments for the ^^Fe+"^*Ti system. 
From the analysis of velocity and energy spectra we 
conclude that no clear evidence of the action of a fis- 
sion barrier could be found. Either fission channels are 
not favoured, or other processes obscured them, like ad- 
ditional evaporation stages or the contribution of many 



mother nuclei rather different in mass. The most relevant 
result is the manifestation of high- velocity tails, which we 
interpreted as possible indications of a preequilibrium ex- 
pansion phase. 



B. Nuclear-model calculations 

We had some hints that very highly excited systems 
are formed even in the ^^Fe-|-p interaction, but we could 
not extract quantitative values directly from the exper- 
iment. We could not recognise the presence of a fission 
barrier, but a more complete analysis is required to ex- 
clude that solely compound-nucleus decays are sufficient 
to explain the Hght- fragment production. Thus, we wish 
to carry out a complete reconstruction of the whole reac- 
tion process and compare the ensemble of experimental 
results with the calculations. 

Henceforth, we will restrict to the ^^Fe-|-p system. In 
particular, we will discuss two possible descriptions for 
the dominant process of light-fragment formation: either 
a series of fission-evaporation decays from a compound 
nucleus, or a fast break-up of a diluted highly excited 
system, in line with a multifragmentation scenario. To 
fix the initial conditions for the two decay models, we 
previously need to calculate mass, charge and excitation- 
energy distributions of hot remnants, as these quantities 
are not observable in the experiment. 



1. Calculation of the excitation energy of the hot collision 
remnants 

The initial non-equilibrium phase of the interaction 
^^'Fe+p was described in the framework of the in- 
tranuclear cascade-exciton model developed by Gudima, 
Mashnik and Toneev [231 • The model describes the in- 
teraction of an hadron or a nucleus traversing a heavy 
ion, considered as a finite open system, composed of 
two degenerate Fermi gases of neutrons and protons in a 
spherical potential well with diffuse surface. The interac- 
tion, pictured as a cascade of quasi-free nucleon-nucleon 
and pion-nucleon collisions, produces high-energy ejec- 
tiles, that leave the system, and low-energy particles that 
are trapped by the nuclear potential. As many holes as 
the number of intranuclear collisions are produced in the 
Fermi gas. The number of trapped particles and the num- 
ber of holes (or excitons, without distinction) determines 
the excitation energy of the so-called "composite nucleus". 

Hot remnants are often treated as equilibrated or par- 
tially equilibrated systems, both in the case of compound- 
nucleus formation and at a freeze-out state. Thus, an 
additional thermaHzation process might be necessary to 
describe the transition from the initial non-equilibrium 
phase of the collision to the equilibrium phase governing 
the decay. Following the hypothesis of the preequiHb- 
rium exciton model, the intranuclear cascade continues 
to develop through the composite nucleus by a sequence 
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of two-body exciton-exciton interactions, until equilib- 
rium is attained. Two kinds of decay characterize the 
composite nucleus: either the transition to a more com- 
plicated exciton state, or the emission of particles into 
the continuum. While in GrifHn's model [^ all decays 
are equiprobable, successive developments proposed more 
elaborate descriptions of the competition between decay 
modes. However, when the conditions of the interaction 
lead to multifragmentation, the evolution of the compos- 
ite nucleus is more compHcated, as the system is sup- 
posed to expand. In the course of the expansion process, 
an intense disordered exchange of charge, mass and en- 
ergy among its constituents is expected. The density of 
nuclear matter evolves to a more dilute state, the freeze- 
out, at which breakup occurs. Sophisticated thermal- 
expansion models were specifically developed to describe 
this thermalization process (lftl21|. Nevertheless, in our 
calculation we were less specific and we adopted Blann's 
preequilibrium exciton model J2a| , independently of the 
deexcitation scenario. 

In fig. El we present a calculation of the hot-fragment 
distribution generated in the intranuclear cascade with 
and without the inclusion of a preequilibrium stage, re- 
spectively. More quantitatively, in fig. El the projec- 
tions of the distribution are shown as a function of the 
excitation-energy-to-mass ratio and mass. We observe 
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FIG. 14: Calculated production of the hot fragments after the 
intranuclear cascade (modelled accordin g to |63 |1 and a pree- 
quilibrium stage (simulated according to |S3l)i ^ot the ^^Fe-|-p 
system. The cross sections of the hot fragments are shown 
as a function of the excitation energy per nucleon E* /A of 
the source (bottom-left) and the mass number (top-right). 
The mass distribution of cross sections of the hot fragments 
is compared to the experimental final-residue production. 



that preequilibrium is particularly effective in evacuating 
part of the excitation energy and widening the distribu- 
tion as a function of the mass. The hot-fragment mass 
distribution is compared to the measured production of 
the final residues, in order to indicate the extension of 
the deexcitation process. 

When preequilibrium is suppressed, the energy per nu- 
cleon available for the deexcitation largely exceeds 2.5 
MeV, a value that corresponds to the temperature of 
around 5 MeV, for a fully thermalised system. In this 
case, the multifragmentation regime is accessible. If pree- 
quilibrium is included, the average excitation of the sys- 
tem extends still right up to the expected threshold for a 
freeze-out state. Hence, according to this previous step 
of the calculation, oscillations in direction to break-up 
decays might be possible. Indeed, at these excitation en- 
ergies break-up channels are expected to be still in com- 
petition with compound-nucleus decay. We will proceed 
to evaluate the extent of this competition by the use of 
deexcitation models. 



FIG. 13: Hot-fra gme nt distributions generated in the in- 
tranuclear cascade ^^, in the case of exclusion (left) and in- 
clusion (right) of a preequilibrium stage |6fil |. for the ^ Fe-|-p 
system. Straight lines define constant values of excitation en- 
ergy per nucleon, indicated in MeV. Five selected isotopes 
correspond to the most probable mass and nuclear charge for 
a given excitation energy per nucleon. 



S. Sequential fission-evaporation decay 

In order to describe the deexcitation process in the 
framework of sequential fission-evaporation decays, we 
applied the code GEMINI 0- Within GEMINI a spe- 
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cial treatment based on the Hauser-Feshbach formalism 
is dedicated to the emission of the Ughtest particles, from 
neutron and proton up to beryllium isotopes. The for- 
mation of heavier nuclei than beryllium is modelled ac- 
cording to the transition-state formalism developed by 
Moretto Q- All asymmetric divisions of the decaying 
compound nuclei are considered in the calculation of 
the probability of successive binary-decay configurations. 
The total-kinetic-energy release in fission originally pa- 
rameterised according to the systematics of Viola |5fll |. 
eq. l|Hl,was replaced by the systematics of Tavares and 
Terranova [^, formulated according to the relation 10, 
and extrapolated for asymmetric splits by the use of the 
conversion (pj. 

We simulated the decay of two possible ensembles of 
hot remnants, those issued directly from the stage of in- 
tranuclear cascade, and those which lost part of excita- 
tion energy and mass in a preequilibrium phase. The 
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FIG. 15: Comparison of the measured mass distributions as 
a function of the mass number for the system ^^Fe+p with 
the results of GEMINI (upper part) and SMM (lower part). 
SMM is more sensitive than GEMINI to the effect of a pree- 
quilibrium phase. 



resulting distributions of final residues are almost indis- 
cernible, revealing that the intermediate-mass fragments 
(especially those around oxygen) are not especially sen- 
sitive to the variation of average excitation energy of 
the system. It might be also pointed out that, when 
very hot fragments are allowed to decay by solely fission- 
evaporation channels, many nucleons and some light clus- 
ters are liberated at the very beginning of the deexci- 
tation, before eventually forming an intermediate-mass 
fragment by fission. When the preequilibrium phase is 
suppressed, this preliminary emission could constitute a 
compensating process. In average, the relation between 
energy loss and mass loss could be similar in the two pro- 
cesses, and lead to analogous results. The difference is 
only conceptual, as the preequilibrium acts on a system 
still evolving toward thermalization, and particle evapo- 
ration is connected to a completely thermalised system. 
Only lithium and beryllium revealed a visible enhance- 
ment in the yields with the increase of average excitation 
energy. 

The result of the model calculation, compared with the 
measured cross sections is presented in the upper side of 
fig-El The evaluation of the heavy-residue cross sections 
is consistent with the experimental data, but a sizeable 
underestimation of the production fails to reproduce the 
intermediate-mass region. Especially the production of 
the residues populating the characteristic hollow in the 
mass distribution reveals to be generally underestimated 
by the calculation. To complete the comparison, we turn 
now back to the first key observable found in our exper- 
imental investigation: the velocity spectra of light frag- 
ments. In the first row of fig.^Jthe experimental spectra 
of ^Li, ^°B, ^^C and ^^C are shown, together with their 
velocity reconstruction (solid line). Within GEMINI, all 
decays are decorrelated in time and when more fragments 
are produced they do not interact in the same Coulomb 
field. Binary compound-nucleus emission is connected 
with a restricted range of heavy sources close to the pro- 
jectile mass, refiected in the small width of the Coulomb 
peaks, as shown in the second row of fig. El This feature 
characterizes only the formation of the lightest fragments 
and disappears with increasing mass of the residues. The 
calculations presented in the second row of fig. 1161 should 
not be compared to the experimental data. The effect of 
the Coulomb repulsion involved in the deexcitation and 
disentangled from the smearing effect of the intra- nuclear 
cascade and preequilibrium emission can be appreciated 
in the third row of fig. El where the reference frame 
has been fixed to the center of mass of the initial sys- 
tem formed at the beginning of the fission-evaporation 
process. In the calculation, the transformation of the 
two Coulomb peaks into one single wide hump occurs for 
lower masses than experimentally observed. The model 
generates one single hump in the longitudinal velocity 
spectra of light fragments when a longer evaporation cas- 
cade is involved, and characterized by mainly alpha and 
nucleon emission. Moreover, the total width of the cal- 
culated spectra is narrower than observed. 
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FIG. 16: First row: Experimental velocity spectra (circles) and reconstructed velocity spectra (solid line). Each spectrum is 
drawn in the reference frame corresponding to the measured average velocity value of the fragment considered. (This frame 
corresponds to the "center of mass" frame of the reaction product drawn in fig. HJ Second row: Calculated velocity spectra 
obtained by GEMINI or SMM following INC and the preequilibrium stage, and from SMM following directly INC. Each 
spectrum is drawn in the reference frame corresponding to the calculated average velocity value of the fragment considered. 
Third row: Velocity recoil introduced by the GEMINI or SMM phase alone (recoils by INC and preequilibrium stages not 
included). All spectra are normalized to the unit. 



3. Fast break-up 

We imputed the underestimation of intermediate-mass 
fragment formation to an incomplete description of the 
most highly excited decaying systems when solely fission- 
evaporation deexcitation was considered. In this respect, 
we turned to the Copenhagen-Moscow statistical multi- 
fragmentation model (SMM) [ll|,[l3|, that is the exten- 
sion of the standard statistical evaporation-fission picture 
toward high excitation energies, treated by adding the 
fast simultaneous disassembly of the system as a pos- 
sible decay channel. The hybrid model of intranuclear 
cascade followed by SMM was already appHed in previ- 
ous studies of proton-induced reactions [67. 68] for the 
description of similar experimental data. In the frame- 
work of SMM, the evaporation from the compound and 
compound-like nuclei is included and, therefore, at low 
excitation energies, if the channels with production of 
compound-like nuclei dominate, SMM gives results simi- 
lar to GEMINI. In particular, the statistical cluster evap- 
oration is treated within the Weisskopf formalism, ex- 
tended to the emission of nuclei (in their ground state 



or available excited states) up to ^*0 J23|. On the other 
hand, when very high excitation energies are reached in 
the collision, the system is assumed to be diluted and 
to have attained the freeze-out density pb- In previous 
studies [l3| Pb was calculated to evolve as a function 
of the excitation energy per nucleon toward an almost 
asymptotic value equal to 1/3 of the ground-state den- 
sity po for high excitation energies {E* /A > 5 MeV). In 
the present calculation, an energy-dependent free volume 
is used to determine the probability for different break- 
up partitions. On the other hand, for the calculation of 
the Coulomb interaction among fragments the freeze-out 
density pb is introduced as a fixed quantity, equal to the 
asymptotic value pb = Po/3. According to the physical 
picture, when the region of phase (spinodal) instability 
is reached, at least partial thermodynamic equilibrium 
is expected and the fragment formation takes place ac- 
cording to chaotic oscillations among different break-up 
configurations, from event to event. In SMM, within the 
total accessible phase space, a microcanonical ensemble 
of all break-up configurations, composed of nucleons and 
excited intermediate-mass fragments governs the disas- 
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FIG. 17: The hatched areas represent portions of the residue 
production calculated with SMM, subdivided according to dif- 
ferent multiplicities of intermediate-mass fragments (having 
A>4). The total production measured experimentally (dots) 
and calculated (solid line) is superimposed for comparison. 
The calculation disregards preequilibrium in the left diagrams 
and includes preequilibrium in the right diagrams. 



sembly of the hot remnant. The probability of different 
channels is proportional to their statistical weight. Sev- 
eral different break-up partitions of the system are pos- 
sible. 

In fig. 1 151 the calculation based on SMM reveals to 
better describe the reaction in comparison to GEM- 
INI. It should also be observed that the production of 
intermediate-mass fragments is sensitive to the excita- 
tion energy of the source. A more detailed view on the 
reconstruction of the reaction mechanism is presented in 
fig. El where the multiplicities involved in the fragment 
formation are investigated. The major cross sections are 
fully determined by evaporation decays. This is true for 
the calculation where preequilibrium is included. On the 
contrary, when preequilibrium is excluded, a depletion of 
the heavy evaporation residues arises as a result of the ex- 
cessive enhancement of higher-multiplicity modes (clus- 
ter emission and multifragmentation) . The intermediate- 
mass fragments are almost totally produced in binary de- 
cays of mainly break-up character. Multifragment emis- 



sion channels (with multiplicity mainly equal to three) 
have a minor contribution in the decay of the thermalised 
system but, when preequilibrium is disregarded, their in- 
cidence is in strong competition with binary splits. We 
might conclude that a consistent description of the pro- 
duction of fragments could be found in between these two 
approaches. 

We can extend the investigation to the excitation en- 
ergies connected to the production of fragments with dif- 
ferent masses. According to the calculation presented 
in fig. El (now performed only including the preequilib- 
rium phase), intermediate-mass fragments are almost all 
formed in the decay of highly excited remnants, with ex- 
citation energy per nucleon above 2.5 MeV. 

An experimental indication of how multiplicity is re- 
lated to the excitation energy is suggested by fig. |HI and 
the upper diagram of fig. El where the ratio of the 
sum of the individual production cross section in the 
two reactions and the total cross-section ratio are com- 
pared. The latter is calculated according to the model of 
Karol [2^. The production cross-section ratio scales with 
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1 MeV < E*/A < 2.5 MeV 



E*/A > 2.5 MeV 
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FIG. 18: Different portions (hatched areas) of the residue 
production calculated with SMM are selected according to 
different ranges in the excitation energy per nucleon E* /A of 
the source. The calculation is performed only for the case of 
inclusion of preequilibrium. The total production measured 
experimentally (dots) and calculated (solid line) is superim- 
posed for comparison. 
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FIG. 19: Upper part. Experimental data on the production 
cross-sections ratio as a function of the mass number for the 
reaction ^'^Fe+"''*Ti versus the reaction ^^Fe+p. The data 
are compared with the ratio of total nuclear cross sections 
for the two reactions, calculated according to the model of 
Karol |2fl|. Lower part. SMM calculation of the probabil- 
ity for the formation of a residue as a function of the mass 
number and the multiplicity. The multiplicity is intended as 
the number of projectile-like residues heavier than an alpha 
particle produced in one collision. 



the total cross-section ratio only in the region of higher 
masses, presumably coming from more peripheral colli- 
sions, while it deviates for lighter masses. The deviation 
must be related to the different mean multiplicities in 
the two reactions: lighter masses are more populated in 
the ^Spg^nat-jj reaction, and great part of this increase 
might be related to higher multiplicities. The observation 
of the gradual increase of multiplicity with the excitation 
of the system, verified in the calculation of different yield 
spectra associated to different energy ranges of hot rem- 
nants, could be followed further when extending to the 
^^Fe+"^*Ti reaction. As pictured in the lower plots of 
fig. El this behaviour is well reproduced by calculations 
with SMM. In this respect, the light-residue production 
characterizing the ^^Fe-|-p system might just be inter- 
preted as the early onset of the process that will govern 
the decay of the sepg^nat-j^ gyg^gjjj_ 

To conclude this section on model calculations, we fo- 
cus once more on the velocity spectra of the light frag- 
ments shown in fig. El We already discussed the diffi- 
culty to combine the wide shapes of the velocity spectra 
and their mean values with a fission barrier. We inferred 
that the extensions of the velocity distributions to very 
high velocities might reflect higher kinetic energies than 
an asymmetric flssion process could release. Consistently 
with this expectation, on the basis of model calculations 
we could connect the production of light residues to very 
high excitation energies of the source. Above around 
2.5 MeV of excitation energy per nucleon, the process of 
light-residue production is still presumably dominated by 



binary decays, but the contribution of disintegration in 
more fragments is not excluded. In this case, parts of the 
distribution corresponding to smaller velocities should be 
more populated than in a purely binary split. In the rep- 
resentation of fig. El this contribution would fill more 
central parts of the spectra when lower asymmetry char- 
acterizes the break-up partition. In a statistics of events 
where three about-equal-size fragments are produced si- 
multaneously, the velocity spectrum of any of them will 
be Gaussian-like. If three fragments are produced, of 
which two are considerably lighter than their heavy part- 
ner, the velocity spectra (Tr of the two light ejectiles will 
be double humped. Another contribution in populating 
lower velocities could be associated to different break- 
up configurations, where the partner or the partners of 
the light residue have different masses. In this case, 
the spectrum is the folding of several binary-like com- 
ponents characterized by different spacing between the 
two maxima and different widths around the maxima, all 
this resulting in the superposition of two (backward and 
forward) triangular-like distributions that could eventu- 
ally merge in a general bell-shape. We might also con- 
sider standard evaporation cooling down the break-up 
residues, emitted in some excited states. In this case, 
the secondary "slow" emission process operates outside 
of the common Coulomb field of the fragmenting rem- 
nant, and would produce a general widening of the spec- 
trum around its maxima. As portrayed in the second 
row of fig. El SMM describes very consistently the ex- 
perimental spectra. In the third row the effect of the 
Coulomb interaction and eventually the expansion is il- 
lustrated by referring to the center of mass of the system 
formed right after the intra-nuclear cascade and the pree- 
quilibrium, if included. SMM calculates the Coulomb in- 
teraction between fragments by placing them inside the 
freeze-out volume ph- Contrarily to GEMINI, it takes 
into account different positions of the fragments, includ- 
ing two-body and many-body partitions. Some multi- 
fragmentation channels may resemble a two-body process 
even at relatively high excitation energy. These channels 
can also include additional small fragments which may 
look like evaporation ones. However, these additional 
small fragments can essentially change the Coulomb in- 
teraction in the volume and the thermal energy in the 
system, and influence the kinetic energies of the main 
two fragments. The binary character characterizing the 
experimental results is properly reproduced and the ve- 
locity distributions calculated with SMM are wider than 
those obtained from GEMINI. As shown in flg. El the 
gradual filling of the center of the spectra could be re- 
lated to different break-up configurations and to possible 
multibody disintegration. The abrupt change of shape in 
passing from ^^C to ^^C (correctly reproduced by SMM) 
might be related to the more favoured evaporation chan- 
nel toward the formation of ^^C, that could collect several 
different decay processes and evaporation decays from 
neighbouring nuclei. This preferential decay toward ^^C 
smears out any binary character of the spectrum. As 



19 



Normalized yield 

cTf(v„)/a" 
(SMM) 




-2 2 
v., [cm/ns] 



FIG. 20: Contribution of different multiplicity channels 
(A>4) to the velocity spectrum of ^Li, as calculated by SMM. 
The representation is the same as in the second column of 

fig.nn 



the break-up configuration varies with the mass and the 
charge of the end-product, it varies also with the excita- 
tion energy available for the disassembly of the hot rem- 
nant. This was evident in fig. El where we compared 
the production cross sections of the residues, and it is 
also evident in fig. El by analysing the velocity spectra. 
It is evident that the suppression of the preequilibrium 
induces a smearing effect on the spectra. This effect is 
dramatic for ^Li as the double-humped spectrum is com- 
pletly smeared out in one large single hump with a fiat 
top, and it becomes similar to the velocity spectrum of 
^Li produced in the even more excited ^^Fe+ "^'Ti sys- 
tem. 



V. CONCLUSION 

The mechanisms of light-particle emission from two 
systems, ^"^Fe-l-p and ^'^Fe+ "'"'Ti at 1 A GeV, have been 
investigated. The latter was regarded as a baseline for 
very high-energy processes (multifragmentation) . We fo- 
cused mainly on the proton-induced reaction. The cer- 
tain understanding that we could attain from the analysis 
of the ^^Fe-l-p system is that light residues are produced 
in the decay of highly excited remnants. Furthermore, 
from a more quantitative discussion, we inferred that the 
emitting source should also be heavy and close to the pro- 
jectile mass. The magnitude of the Coulomb repulsion, 
together with the very high formation yields even sug- 
gested that an asymmetric break-up process, hardly con- 
nected to asymmetric fission or statistical cluster emis- 
sion, might be the favoured channel of light-residue pro- 
duction. These findings were derived from experimen- 
tal observables like the isotopic cross sections measured 
for the whole ensemble of the residues, and the veloc- 
ity distributions of the emitted fragments in the projec- 
tile frame along the beam axis. Especially the shape of 
the velocity spectra offered us a microscopic insight into 
the mechanisms of light-particle emission. In analysing 
the features of the velocity spectra, we failed in describ- 
ing the kinematics within a general systematics of fis- 
sion total-kinetic-energy release. A complete simulation 
of the whole reaction process, where sequential fission- 
evaporation decays govern the deexcitation, could not 
consistently describe the gross experimental features of 



the decay. 

We suggested that the characteristics of the kinematics 
and the production of light particles could carry indica- 
tions of fast asymmetric splits. A novel description of 
the complete reaction process, including channels of fast 
break-up decays revealed to be more adapted in depicting 
the decay of the most highly excited remnants, and was 
compatible with the high yields for light residues and the 
complex shapes of the velocity spectra. Encouraged by 
this consistency and, first of all, on the basis of previous 
theoretical and experimental results (see references in the 
section ITvll . we suggested that protons at incident ener- 
gies of 1 A GeV traversing iron-like nuclei can introduce 
very high thermal excitation energy per nucleon in the 
system, even above 2.5 MeV. Such a thermal excitation 
could lead to attain freeze-out conditions. 

We assume that if the excitation energy is high but 
not sufficient to lead to freeze-out, the system might still 
expand, but without reaching the break-up phase, the ex- 
pansion subsequently reverses to compression in a path 
toward the formation of a compound nucleus . This be- 
haviour might occur before the nucleus cools down by 
fission-evaporation decays. On the contrary, when the ex- 
citation energy is just sufficient to access break-up chan- 
nels, partitions with low multiplicity of intermediate- 
mass fragments and high asymmetry are favoured: the 
decay results mainly in the simultaneous formation of one 
heavy residue, with mass close to the hot remnant, and 
one or more light clusters and nucleons. As an extreme 
case, two fragments rather asymmetric in mass may be 
formed in the same fast break-up process. The forma- 
tion of light fragments in the ^^Fe+p reaction could be 
explained by this picture. 



VI. 
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APPENDIX A: THE VELOCITY 
RECONSTRUCTION 

1. Demonstration of the relation l(3JI 

We describe hereafter the derivation of the equa- 
tion ||SJ, which connects the measured spectra 
dr(w|| )/dw|| as a function of the longitudinal velocity com- 
ponent v\\ in the center-of-mass frame to the cross-section 
variation in velocity space in the center-of-mass frame. In 
a general case, the latter distribution is not isotropic, but 
a function of the absolute velocity v, the polar angle from 
the beam direction 0, and the azimuthal angle around the 
beam axis ip. It will be denoted as d^a/{dvd^), where 
fl is the solid angle. The velocity component orthogonal 
to the beam axis is v±. The contribution to the experi- 
mental yield in the interval [w||,W|| + Aw||] is obtained by 
integrating v± in the slab orthogonal to the beam axis : 



cg(^ll) 
dv\\ 



——- v± dv± dip 
dv 

1 dV 



v'^ dv dil 



vj_ dv± dip 



(Al) 

(A2) 



2. Reversibility of the relation ({S)) 

The physical quantity of interest in equation ^ is the 
term da/dv, describing the variation of the cross sec- 
tion a(v) as a function of the absolute velocity v in the 
center-of-mass frame, whereas the measured quantity is 
the left-hand term dX(w|| )/dw|| , representing the variation 
of the apparent cross section as a function of the longitu- 
dinal velocity component v\\ in the center-of-mass frame. 
In principle, equation (O could not be inverted in an un- 
ambiguous way for general shapes of the dcr/du function. 
However, for the restricted shapes describing the data, 
this inversion becomes possible. This is particularly the 
case if this function is supposed to decrease monotoni- 
cally to at large v and if dl(t'||)/di;|| also follows the 
same behaviour at large \v\\\, as it is evident from fig. El 

The inversion procedure can be described as follows. 
Let's consider a given bin in longitudinal velocity defined 
by the interval ["yitiWn + Afy]. The yield for this bin is 
(dJ(i'||)/dw||)/Aw||, while the corresponding integral over 

V in equation ^ extends from v\\ to 



^11 



This 



domain is depicted by the thick segment in fig. [2] Let us 
assume that the values of the function dcr/dv are known 
over this interval and that they comply to equation JSJ . 



^"f" 



V (V||-AV|| f+vl 



V^ 



+vf 



For the orthogonal velocity integration the lower limit is 
and the higher limit is related to the angular accep- 
tance of the spectrometer. Since the latter is not neces- 
sarily circular, it can depend on p and will be denoted as 
a{p). The maximal orthogonal velocity may be derived 
from the Lorentz transformation of the momentum and 
it reads : v±{ip) = ^(u -\- v\\)a{ip), where u and 7 are 
the velocity and the Lorentz factor of the center of mass 
in the laboratory frame, respectively. Introducing these 
limits in the integration, we write : 



dvw 



-iv) 



1 dV 



v^ dv dQ 



vj_ dvj 



dp 



(A3) 



Changing the integration variable from v± to v 



Vu + v^ we obtain 



dz;ii 



/vf.+vl{v) 



1 d^cr 



V dvdfl 



dv 



l^-iil 



dp 



(A4) 



In the case of an isotropic velocity distribution, 
d^a/{dvdn) reduces to {1/ An) (da/dv) which, substi- 
tuted in eq. (|A4|I gives eq. (EJ. 



FIG. 21: Integration domains of eq. Q. 



We consider now the W|[ bin located between U|| — Ai)|| 
and I'll , for which the integration extends from vu — Av|| 

and w(w|| — Aw||)2 + vj^ as shown by the thin segment 

above the axis in fig. |21 It can be seen that this seg- 
ment has a large overlap with the previous one, provided 
Aw|| is small enough. If this were not the case, the inver- 
sion procedure would even be simplified as da/dv would 
be directly proportional to the yield divided by the in- 
terval length, provided that the yield has a low variation 
over Aw|| . In the case of overlap of the two integration 
segments, as shown in the figure, the variation of the 
yield comes only from the values of da/dv at the edges. 
If the value is known on the right non-overlapping ex- 
tremity, the variation between the two adjacent z;|| bins 
delivers the value on the left non-overlapping extremity. 
The procedure can be continued for lower v^\ bins, fixing 
the values of the function for decreasing values of v. 

So far, no specific assumption has been made except 
that the function da/dv is known over a given inter- 
val. This can be practically achieved if one assumes that 
da/dv vanishes at large v values and that, as a conse- 
quence, also the yield drops. By considering a large v^^ 
value for which the yield is null, we can take a null da/dv 
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over the corresponding interval and start the procedure 
of reversion. This prescription for the starting point can 
also be extended to regions where the yield does not fade : 
for W|| ;2> v± the length of the integration interval de- 
creases as v]_/{2v\\), which becomes small compared to 
the characteristic variation length oida/dv. In this case, 
the latter can be assumed constant over the interval and 
its value deduced straightforward from eq. ^. The de- 
pendence of v± on ifi only slighly changes the procedure, 
while the scheme remains the same. 

The yields measured for the forward emission (w|| > 0) 
are expected to differ from those associated to the back- 
ward {v\\ < 0) emission. Nevertheless, in the ideal case 
of a perfectly isotropic emission with respect to the cen- 
ter of mass, the resulting cross sections cr(u|| > 0) and 
ct{v\\ < 0) restricted to only-forward and only-backward 
emission, respectively, should be identical. The differ- 
ence |o'(w|| > 0) — cr(w|| < 0)1 can be an indication of 
the uncertainty introduced in the extraction of the cross 
section a{v) by the assumption of isotropic emission. 



APPENDIX B: ISOTOPIC CROSS SECTIONS 



In the subsection III Bl we described how cross sections 
could be extracted from the measured longitudinal ve- 
locity spectra. The results are presented in tabled with 
the statistical uncertainties. In appendix El the numeri- 
cal tool used in our analysis was presented. It should be 
observed that our velocity-reconstruction method allows 
to obtain cross sections for isotopes of which at least a 
half of the longitudinal velocity spectrum is measured. 
In this case, some needed parameters like the mean re- 
coil velocity or the width of the distribution could be 
extrapolated from neighbouring isotopes. On the other 
hand, the whole procedure is valid up to a certain ex- 
tent: due to the assumption of isotropic emission, ideal 
cases should result into equal (equal area and equal cen- 
troid) absolute-velocity distributions deduced from the 
forward and the backward part of the measured longitu- 
dinal velocity distributions. Deviations from this ideal 
case derive either from the physics of the reaction pro- 
cess, that could differ from a purely isotropic emission , 
or by the lack of statistics in some parts of the spectrum, 
resulting in complicating the convergence of the numer- 
ical calculation. This leads to results that fluctuate by 
10% in the average. We take this value as the statistical 
uncertainty (and not simply the statistics of counts) . 



TABLE I: Spallation and fragmentation residue isotopic cross 
sections measured in this work for the formation of Li, Be, B, 
C, N and O in the reaction ^^Fe+p and •■^'^Fe+"''*Ti, respec- 
tively. 



Isotope 


^''Fe+p, a [mb] 


^^Fe+""*Ti, a [mb] 


^Li 


14.89 ± 1.5 


128.50 ± 12.9 


^Li 


3.06 ± 0.3 


103.46 ± 10.3 


^Be 


3.09 ± 0.3 


62.28 ± 6.2 


^Be 


2.11 ± 0.2 


29.88 ± 3.0 


10 g 


2.03 ± 0.2 


35.06 ± 3.5 


"B 


3.82 ± 0.4 


61.58 ± 6.2 


12b 


0.36 ± 0.04 


9.56 ± 1.0 


"C 


1.12 ± 0.1 


16.02 ± 1.6 


12^ 


4.69 ± 0.4 


62.50 ± 6.3 


13 c 


2.76 ± 0.3 


37.59 ± 3.6 


"c 


1.87 ± 0.2 


14.74 ± 1.5 


"N 


0.14 ± 0.01 




14n 


1.25 ± 0.1 


18.67 ± 1.9 


15N 


2.97 ± 0.3 


37.33 ± 3.7 


16n 


0.33 ± 0.03 


5.04 ± 0.5 


i^N 


0.11 ± 0.01 


2.44 ± 0.2 


18N 




0.01 ± 0.001 


140 


0.01 ± 0.001 




15q 


0.33 ± 0.03 


5.91 ± 0.6 


16(3 


2.78 ± 0.3 


36.62 ± 3.7 


"0 


1.46 ± 0.1 


15.58 ± 1.5 


18 


0.75 ± 0.08 


9.12 ± 0.9 


19 


0.13 ± 0.01 


1.99 ± 0.2 



The systematic uncertainties are in general very small 
in FRS measurements of spallation residues. Indeed they 
rise to considerably high values when the measurement is 
dedicated to fragments having very high velocities in the 
projectile frame. This is the case of very light fragments 
emitted in fission-like events or in break-up processes. 
The largest source of uncertainty is the angular accep- 
tance. Heavy residues, close to the projectile mass are 
emitted very forward, and the angular acceptance is close 
to 100%. On the contrary, light fragments are strongly 
affected. The multiplicity of the intermediate-mass frag- 
ments could not directly be measured. From physical 
arguments we could safely infer that light fragments are 
emitted in events with multiplicity (of fragments with 
A>4) prevalently equal to two. Indeed we could not 
exclude the possible contribution of higher multiplicity 
processes. We estimated the systematic uncertainty to 
be up to 30%. 



It might be remarked that the greatest contribution 
to the total uncertainty comes from the systematic un- 
certainty. On the other hand, the cross-section ratios of 
different nuclides are very consistent as they are related 
to small statistical uncertainties. 
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